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Abstract 

The discrete heat equation is worked out in order to illustrate the search of sym- 
metries of difference equations. It is paid an special attention to the Lie structure 
of these symmetries, as well as to their dependence on the derivative discretization. 
The case of g-symmetries for discrete equations in a g-lattice is also considered. 



1 Introduction 

As it is well known Lie point symmetries were introduced by Lie for solving differential 
equations, providing one of the most efficient methods for obtaining exact analytical 
solutions of partial differential equations j0J. The interest for discrete systems in the last 
years has led to extend the Lie method to the case of discrete equations @-|§]. 

A general difference equation, involving one scalar function u(x) of p independent 
variables ) evaluated at a finite number of points on a lattice will be 

written in the form 

E(x, T a u(x), T hi A x .u(x), V A.,,A,,„(.r). . . .) = 0, (1.1) 

where the shift operators T a ,T b \ T Clj are defined by 

T a u(x):={T^T^---T^u(x)} ni _ , a = ( 0l> a 2 , . . . , a p ), i = l,2,...,p, 
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with a,i, mi, rii, (rrii < rii), fixed integers, 

and <7j is the positive lattice spacing in the uniform lattice of the variable X{ (i = 1, . . . ,p). 
The other shift operators T ' , T CiJ are defined in a similar way. The difference operators 
A x . are defined so that in the continuous limit turn into partial derivatives. 

In the following we will make use of the approach presented in ||, based on the 
formalism of continuous evolutionary vector fields . The infinitesimal symmetry vectors 
in evolutionary form for the difference equation of order N given in ( |1 . 1| ) take the general 
expression 

X e = Qdu = (^2 6 T a u, a x , a t )T b A Xi u - <j>(x, T c u, a x , a t )^j du, (1.2) 

where £,i(x,T a u,cr x ,a t ) and (p(x,T c u,a x ,a t ) are operator valued functions which in the 
continuous limit become the functions u) and 0(x, u), respectively, giving rise to Lie 
point symmetries. 

The vector fields X e generate the symmetry algebra of the discrete equation (|1 . 1|) , 
whose elements transform solutions u(x) of the equation into solutions u(x). The N— th 
prolongation of X e must verify the invariance condition 

pr N X e E\ E=0 = 0. (1.3) 

The prolongation pr N X e is 

pr N X e = J2 T a Qd T a u + J2T bl Q x *d TblAxU + £ T^Q^d^ A , u + . . . (1.4) 

where summations in ( |1.4j ) are over all the sites present in ( |1.1| ). The symbols Q Xi , Q XiX \ . . . 
are total variations of Q, i.e., Q Xi = A^Q, Q XiX i = A Xi A x .Q, ■ ■ ■ , defined by 

Alf(x,u(x),A x u(x), ■■■) = -[f{x+a,u(x+a), (A x u)(x+a), . . .)-f(x,u(x),A x u(x), . . .)], 

a 

while the partial variation A x . is 
A x f(x,u(x), A x u(x), . . .)=^[f(x+a,u(x), (A x u)(x), . . .)-f(x,u(x),A x u(x), . . .)]. (1.5) 

The solutions of ( |1.3|) gives the symmetries of equation ( |1 . 1|) when using the difference 
operator ( |1.5| ). The determining equations for & and are obtained by considering linearly 
independent the expressions in the discrete derivatives T a A x .u, T b A x . Xj u, .... The Lie 
commutators of the vector fields X e are obtained by commuting their first prolongations 
and projecting onto the symmetry algebra Q . 

Since we will restrict here to linear equations we can assume that the evolutionary 
vectors (|1.2j ) have the form X e = (Xu)d u , where 

X = £ T a , a x , a t )A x . - <f>(x, T a , a x , a t ). 
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The operators X, in general, may span only a subalgebra of the whole Lie symmetry 
algebra ||. 

In Ref. H the symmetries were obtained using the above mentioned difference operator 

A Xi ee A+ = ^± (1.6) 

which, when crj — > 0, goes into the standard right derivative with respect to Xj. Since 
other definitions of the difference operator can be introduced 0, one would like to show 
that the algebraic structure of the symmetries is independent on the choice undertaken. 
In this work we will also use the left derivative 

1 - T- 1 

A" = (1.7) 

and the symmetric derivative (which goes into the derivative with respect to xi up to 
terms of order of) 

rp _ rp-l 

Af.. = — (1.8) 



2tr, 



In the following we will see that by an appropriate definition of the Leibniz rule we can 
construct Lie symmetries, in principle, for any difference operator. In Section Q we will 
introduce this procedure on the example of the discrete heat equation. We shall study 
separately the cases of the discrete derivatives (|1.6|) - (|1.8|) , as well as that of a g-derivative. 
In this way we will get different representations of the same Lie algebra. We conclude 
with some remarks and comments. 



2 Discrete heat equation 

Let us consider the second order difference equation 

(A t - A xx )u(x) = 0, 

as a discretization of the heat equation. Since it is linear, we can consider an evolutionary 
vector field of the form 

X e = Qd u = (rA t + £A x u + fu)d u , (2.1) 

where r, £ and / are (operator valued) functions of x, t, T x T t , a x and a t . The determining 
equation is 

AjQ - Al x Q\ Axxll=AtU = 0, 

whose explicit expression is 
At^A^+A^rA^+A^/^-tA^^A^+A^^A^+A^^/M)]^^^ = 0. (2.2) 
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Only when expression ( |2.2p is developed one needs to apply a Leibniz rule and, hence, 
the results will depend from the definition of the corresponding discrete derivative. We 
propose a Leibniz rule having the form 

A, (f(x)g(x)) = f(x)A x g(x) + D x (f(x))g(x), (2.3) 

where D x (f(x)) = [A x , f(x)] is a function of x,T x and a x (similarly for D t (f(t))). 

Using the general rule ( [2.3| ) for an arbitrary discrete derivative we obtain from ( |2.2| ). 



equating to zero the coefficients of A xt u, A t u, A x u and u, respectively, the following set 
of determining equations 

D x (t) = 0, 

DAt) - 2Djn = 0, 

(2.4) 

D t (O-D xx (O-2D x (f) = 0, 
D t (f) - D xx (f) = 

where D xx (f) = D x (D x (f)). Next, starting from ( |2.4| ) we will study separately the cases 
for A ± and A s . 



2.1 Symmetries for right (left) discrete derivatives 

Choosing as in Ref. || || the derivative A + and, consequently, the Leibniz rule 

A + (fg) = fA + g + A + (f)Tg, 

we get from (|2]4]), 

A+r = 0, 

(A+r)T t - 2(A+$T X = 0, 

(Att)T t - {AU)T 2 X - 2(A+ f)T x = 0, 

{Atm - (a+ m = o- 

The solution of (|2.5| ) gives 

t = tP>T 2 + tn + r , 

£ = |x(ri + 2tr 2 )T t T x 1 + i& + £„, (2.6) 
/ = ix( 2 )r 2 T t 2 T- 2 + |tr 2 T t + i^T" 1 + 7 , 

where r , r l5 r 2 , £o> £i an d 7 are arbitrary functions of T x , T t , and the spacings a x and 
cr t . The notation x^ n \ t^ 71 ' is for Pochhammer symbols; for instance 

x^ = x(x — a x ) . . . (x — [n — l)^). 
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By a suitable choice of the functions 7$, and 7 we get the following symmetries || 

p = (A t «)a ttJ 

A = (A x u)d„, 
W = ud u , 

B = (2tTf 1 A x u + xT^ujdu, (2.7) 

D = {2tT t ~ l A t u + xT~ l A x u + i«)<9 U) 

# = (t 2 Tf 2 A t u - a t tTf 2 A t u + txTf 1 ^ 1 A x u 
+ \x 2 T~ 2 u - \a x xT~ 2 u + \tTr l u)d u . 

Let us note that the above discrete symmetries have a well defined limit when a x , at — > 0, 
which leads to the symmetries of the continuous heat equation. Also, it can be checked 
that with this choice, the symmetries fl2.7p close a 6-dimensional Lie algebra isomorphic 
to the symmetry algebra of the continuous heat equation valid for any value of a x ,a t . 

A second choice for the discrete derivative is A~. The Leibniz rule becomes 

A-(fg) = fA-g + A-(f)T- 1 g. (2.8) 
It gives the same results ( p.6[) and ( p.7|) provided we make the substitution T — > T^ 1 . 

2.2 Symmetries for symmetric discrete derivatives 



Next, let us consider the case of the symmetric derivative A s Q1.8Q . The commutator 

can always rewritten by introducing a function f3 x = (3 S {T X ) = 2(T X + T" 1 ) -1 as 

[A x , x(3 x ] = 1. 

This fact will help us in the computation of general commutators, so 

[A x , f(.r)C = (A x f(x)) T x i; + (V;,. -fir) - fir)) (3 X A S X . (2.9) 
From relation ( |2.9| ) we get the explicit Leibniz rule 

A: (f(x)g(x)) = f(x)A x g(x) + [ i- ((T- 1 - l)/(x)) (T x - (/J*)" 1 ) + (A s J(x)) T x \ g(x). 

This formula allows us to write explicitly the determining equations (|2.4j) . Their solution 
is given by 

f = |x (2tr 2 + n + crtTf 1 ^) {(3t)~ 1 f3 s x + t& + 
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where r 2 , r±, t , £1, £o and /o are arbitrary functions of T x , T t , a x and a t . Now, from 
these solutions and ( |2.1j ) we obtain, with a suitable choice of r 2 , tl, r , £1, £o and fo, the 
following symmetries 

P$ = (Atu)d u , 
PI = (A x u)d u , 
W s = ud u , 

B s = (2*# A> + x(3 s x u)d u , (2.10) 
D s = (2t[3 s t A s t u + x(3 s x A s x u + \u)d u , 

K s = ((t 2 (AS) 2 _ tat 2 m 3 A s )A s u + tx( 3f^A x U 

-\xa x \(3 s x fA%u + \x\fc) 2 u + \mu)d u . 

These symmetries close the same 6-dimensional Lie algebra generated by the operators 
(|2.7|), and have a well defined continuous limit. 



2.3 Symmetries for q— derivatives 

In the following we shall extend the preceding method to the case of ^-derivatives and 
g-symmetries. We deal briefly with a g-discretized heat equation where the g-difference 
operator is defined by 

a«= , * (r,-i), 

[q x -l)x 

with the help of a g-shift operator 

T x = e«* x9 *, T x f(x) = f(q x x). 
In this case we have the commutator 

[Al,x}=T x . 
If we look for a function (5 X {T X ) satisfying 

[Al,(3 x x] = l, 

a formal solution is 

Px{T x ) = (q x - 1) 



T x - 1 

Thus, we can perform a change of basic operators 

{x,T x } ^ {x,A q x }, x = {3 x (T x )x, 
so that, formally, we can express any function f(x, T x ) as 

f{x,T x ) = f(x,Al). 
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In this way the determining equations Q2.4| ) take the form 

r x = 0, 

£t ~~ £,xx — 2 fx — 0, 
ft ~ fxx — 0. 

Therefore, we can give a set of solutions that have a similar appearance to the classical 
symmetries 

P^(A q u)d u , 
P? = (Alu)d u , 

w q = ud u , 

B q = (2(3 t tAlu + f3 x xu)d u , 
D q = (2A*A| it + ^xA' u + ^ 

K q = {-itt 2 A q t u + (3 x (3 t txA q x u + ^ lx x 2 u + iftt m)<9 m , 

where 7 X = ^ 9 z}^ x f x (g- 1 )(- 1 +^^) ff ox replace x by £) . 

L J- x t ] L q x l x — 1 J 

3 Conclusions 

We remark that the key point in obtaining the explicit determining equations ( |2.4| ) and, 
consequently, the discrete symmetries is the use of the Leibniz rule defined in ( [2.3|) . Of 
course, this approach is not the only possibility; in fact, it must be checked whether it 
works correctly or not for each case under study. The choice of a commutator in order to 
define the Leibniz rule implicitly leads us to Lie symmetries, since the natural algebraic 
structure will be given also in terms of commutators. 

Some of the above result deserves some comments. The symmetries associated to the 
symmetric derivatives (Section |2.2|) include functions f3f of T t (T x ) that can only be 
understood as infinite series expansions. Therefore, not all the symmetries ( 2.10|) have 



a local character, in the sense that they are not (finite) polynomials in the operators 
T^ 1 , T^ 1 . Note that although these discrete symmetries give rise to the classical sym- 
metries in the limit a x — > 0, o t — > 0, one of them, K s , also includes surprisingly a term in 
(A t ) 2 , which vanishes in the continuous limit since it contains a factor of. 

Something similar happens with the g-symmetries (Section |2.3| ): they have also a 
highly non-local character. The origin of this unpleasant feature is that the basic com- 
mutator [A|,^x] = 1 needs a non-local function f3 x {T x ). If we want to investigate 
local symmetries it is necessary to use another commutator free of this problem. For 
instance, we could take as starting point the (/-commutator: [A q x , xT~ 1 ] gx = 1, where 
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[A, B} qx = AB — q x 1 BA. However, this treatment would lead us to a g-algebra |S| which 
is out of our present scope. 

Let us insist that the procedure here exposed can be straightforwardly applied to other 
discretizations such as the wave equation [[UJ or even equations including a potential 



term as long as we keep inside the field of linear equations. Non-linear equations needs 
additional improvements in order to have reasonable determining equations. 
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